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Abstract 

Electrochemical etching of semiconductors, apart from many 
technical applications, provides an interesting experimen- 
tal setup for self-organized structure formation capable 
e.g. of regular, diameter-modulated, and branching pores. 
The underlying dynamical processes governing current trans- 
fer and structure formation are described by the Current- 
Burst-Model: all dissolution processes are assumed to oc- 
cur inhomogeneously in time and space as a Current Burst 
(CB); the properties and interactions between CB's are 
described by a number of material- and chemistry- de- 
pendent ingredients, like passivation and aging of surfaces 
in different crystallographic orientations, giving a qualita- 
tive understanding of resulting pore morphologies. These 
morphologies cannot be influenced only by the current, by 
chemical, material and other etching conditions, but also 
by an open-loop control, triggering the time scale given 
by the oxide dissolution time. With this method, under 
conditions where only branching pores occur, the addi- 
tional signal hinders side pore formation resulting in reg- 
ular pores with modulated diameter. 
Silicon monocrystals are the most perfect mate- 
rial mankind has ever created since they are es- 
sential for high integrated computer circuits. Due 
to this high degree of perfection, i. e. corrosion is 
not defect-driven, and "ideal corrosion" is possi- 
ble: Silicon and other semiconductors can be made 
porous by electrochemical etching giving an out- 
standing variety of pore sizes, from nanopores to 
macropores, and geometries, including disordered 
nanoporous, dendritic-like sidebranching pores, and 
pores with modulated diameter. Abstracting from 
the details of the underlying electrochemical pro- 
cesses, the Current Burst Model together with the 
Aging Concept, and accounting the interactions 
between Current Bursts, the generation of differ- 
ent pore geometries, oscillations and synchroniza- 
tion phenomena can be explained including the 
percolation transition to global oscillations. Based 
on the time scale derived from the Aging Concept, 
an open-loop control can be applied to suppress 
sidebranching of pores in a technologically relevant 
regime. 



1 Introduction 

The solid - liquid junction of Silicon and HF - contain- 
ing liquids exhibits a number of peculiar features, e.g. a 
very low density of surface states, i.e. an extremely well 
"passivated" interface [1]. If the junction is biased, the 
IV - characteristics (Fig. 1) in diluted HF is quite compli- 
cated and exhibits two current peaks and strong current- 
or voltage oscillations at large current densities (for re- 
views see [2, 3]). These oscillations have been decribed 
quantitatively by the Current-Burst-Model [4, 5, 6, 7]. 
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Figure 1: The IV- characteristics of the silicon- 
hydrofluoric acid contact shows different phenomena from 
generation of a porous silicon layer (PSL), oxidation and 
electropolishing (OX) and electrochemical oscillations at 
higher anodic bias. 

Perhaps the most outstanding features are the many 
different kinds of pores - nanopores, mesopores, macrop- 
ores, and so on - that form under a wide range of conditions 
in many HF containing electrolytes, including organic sub- 
stances [8, 9]. Despite of an intensive research triggered by 
the finding that nanoporous Si shows strong luminescence 
[10, 11], neither the intricacies of the IV - characteristics 
nor the processes responsible for the formation of pores, 



including their rather peculiar dependence on the crystal 
orientation, are well understood. 

Replacing Si by III/V-compounds, a variety of different 
pore morphologies can be etched; due to different proper- 
ties of the A- and B- surfaces one finds e.g. tetrahedron- 
shaped pores instead of octahedron-shaped pores. For an 
overview over recent results see [12]. But again most of 
the phenomenon can be well understood within the frame- 
work of the Current-Burst-Model which seems to reflect a 
number of quite general properties of semiconductor elec- 
trochemistry. 

2 Experimental Setup and Basic IV 
Characteristics 
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Figure 2: The experimental set-up used for electrochemi- 
cal anodization of semiconductors. 

The basic setup is shown in Fig. 2. Using a four elec- 
trode arrangement a potentiostat/galvanostat is contact- 
ing the sample and the electrolyte, allowing for a well 
defind potential resp. current for the electrochemical dis- 
solution reactions. Since the potentiostat/galvanostat as 
well as the temperature are PC controlled, all relevant 
etching parameters can be controlled in detail. While the 
principal setup remains the same in all experiments, back- 
side contact, front- and/or backside illumination and elec- 
trolyte pumping can be varied as well as cell size (from 
under 0.3 cm up to wafers of 6 in) and semiconductor 
material (Si, InP, GaAs, GaP) including various doping 
levels and crystallographic orientations. In addition, the 
electrolytes (e.g. HF, HC1, H2SO4) and their concentra- 
tions and temperature can be varied. 

3 The Current-Burst Model 

The Current Burst Model [4, 5, 6, 13] states that the dis- 
solution mainly takes place on small spots in short events, 
starting with a direct Si-dissolution, and possibly followed 
by an oxidizing reaction (see Fig. 3). 



Charge 

I 

mean 

Charge\^ _l_ 



> Time 



direct ' ^ ^ 

Si-Dissolution oxidizing Oxide H-Termination Nucleation 2 



Reaction Dissolution 



b) 



Charge 



Charge ^ 



Charge 



J Poretip 



direct 
Si-Dissolution 



oxidizing Oxide H-Termination Nucleation 2 
Reaction Dissolution 



-> Time 



H-Termination 



c) 



d) 



:::: 



Figure 3: Applying a global current density smaller than 
the average current density in a current burst will either 
require very long time constants to keep the mean charge 
density passed in a burst small and to cover the surface 
completely with current lines (a), a statistical arrange- 
ment of current lines with the optimized (smaller) time 
constant of the system, leaving parts of the surface with- 
out current and creating nanopores (b, c), or induces a 
phase separation by rearranging the current lines in areas 
corresponding to macropores (d). 

After these two short processes, the oxide hump under- 
goes dissolution, a time-consuming process which ensures 
at the location of the current burst a dead-time of fixed 
length during which no new burst can start. However, im- 
mediately after dissolution the Si surface has the highest 
reactivity, resulting in a maximal probability of annother 
current burst. Due to H-termination the surface becomes 
passivated, and the probability for bursts decays until it 
reaches the properties of a completely passivated Si sur- 
face, comparable to the situation before the nucleation of 
the first Current Burst. 

This approach allows to connect the average current 
density j with a series of charge and time consuming pro- 
cesses 



En' 



(l) 



For most of these processes the dependence on the choice 
of the etching conditions is known; so by "designing" the 
electrochemical etching conditions (electrolyte, current den- 
sity, applied voltage, temperature, illumination,...) the 
properties and interactions of the basic etching process 
can be tailored to create various pore morphologies. In 
addition (1) directly couples the local average current den- 
sity j to an intrinsic time constant r^. The influence 
of both, current density and intrinsic time constant, can 
be measured independently in experiments. Independent 



measurements can be described with the same set of pa- 
rameters. Within the current burst model pore formation 
from the nm to the fim range as illustrated by Fig. 3 c+d 
is very easy to describe because the separation of the dis- 
solution processes and the surface passivation processes is 
already an immanent property of each dissolution cycle, 
i.e. current burst. 



Figure 4: Schematic view of a "ionic breakthrough" . Due 
to the oxidizing current a roughly semi-spherical oxide in- 
clusion is growing from the tip of channel 1. This is the 
situation of extremely high current (electrochemical oscil- 
lations). 

4 Smoothing effect of Oxide Lay- 
ers and global current oscillations 

Depending on the regime in the IV-curve, Current Bursts 
and their interaction play different roles: For extremely 
high currents, in the oscillation regime, one has a per- 
manent oxide coverage, and due to the high forcing of the 
system current bursts are started at all locations where the 
oxide layer is sufficiently thin for a breakthrough (Fig. 4). 
In this regime, a detailed Monte Carlo study [5, 6, 7] 
including the lateral interaction between Current Bursts 
(overlap of the oxide humps of neighboring Current Bursts) 
has shown that the Current Burst Model quantitatively 
can explain the experimental observations of globally os- 
cillating etching current. Due to a phase synchronization 
of neighboring Current Bursts oxide domains are formed. 
The size of this domains increases with increasing oxide 
generation/reduced oxide dissolution. At a percolation 
point only one oxide domain exists on the sample surface 
with a synchronized cycle of oxide growth and dissolution, 



resulting in a macroscopic oscillation of the external cur- 
rent. The size of the domains as well as the oscillation 
time can be controlled by the chemical parameters. Even 
without global oscillations for all regions of the IV-curve 
where oxide is formed one finds domains of synchronized 
oxide growth which define the length scale for the rough- 
ness of the electrochemically polished sample surface and 
thus lead to a smoothing of the surface. 

5 Passivation Effects: The Aging 
Concept 

While at high current densities the semiconductor surface 
is completely covered with oxide at low current densities, 
most of the semiconductor surface will be in direct contact 
to the electrolyte. It is well known [1] that after chemical 
dissolution the free surface is passivated, i.e. the density of 
surface states reduces as a function of time which increases 
the stability of the surface against further electrochemical 
attack. Schematically the perfection of the surface passi- 
vation and the resulting reduction of the probability for a 
chemical attack as a function of time are plotted in Fig. 
5. For the example of silicon the speed and the perfec- 
tion of passivation of the (111) crystallographic surface is 
larger than for the (100) surface. This selective aging of 
surfaces leads to a self amplifying dissolution of (100) sur- 
faces (which will become pore tips) and a preferential pas- 
sivation of (111) surfaces (which will become pore walls). 
Under optimized chemical conditions with an extremely 
large passivation difference between (111) and (100) sur- 
faces a self organized growth of octahedral cavities occurs 
as is schematically drawn in Fig. 6. The octahedra con- 
sists of (111) pore walls. As soon as the complete surface 
of the octahedra reaches a critical value, it is easier to 
start a new cavity at a (100) tip of the old cavity, since 
the current density in the new, small cavity is larger and 
no surface passivation will occur until the surface again 
becomes to large. This growing mechanism leads to an 
oscillation of both the current through each pore and the 
diameter of each pore as a function of time which is also 
plotted in Fig. 6. As in the case of oxide dissolution an 
internal time constant is related to the pore growth. 

H-Termination H-Termination 




Figure 5: Time constant and perfection of H-termination 
differ strongly for (100) and (111) Si surface orientation. 
This leads to an increased probability of local current nu- 
cleation on (100) surfaces. 
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Figure 6: A consequence of the crystallographic depen- 
dence of the surface passivation in Fig. 5 is the growth of 
octahedra like pores. This pore growth is associated with 
an oscillation of the current through this pore, even for a 
fixed externally applied potential. 

If a global current density smaller than the average cur- 
rent density in a current burst is applied, the aging concept 
becomes essential to explain the generation of macropores: 
At passivated surfaces, the nucleation probability of new 
current bursts is lower than at sites where a current burst 
has taken place before. Therefore the situation shown in 
Fig. 3(c) is much more unlikely than the situation shown 
in Fig. 3(d), so that the whole surface separates in two 
phases with current-carrying pores and passivated surface 
without contributions to the total current. 

6 Consequences on Pore Geome- 
tries 

As described above, two synergetic effects occur as a con- 
sequence of the interaction of Current Bursts: 

• A surface smoothing due to the dissolution of an 
isotropic oxide layer. 

• A strongly anisotropic dissolution of the semicon- 
ductor due to the crystallography dependent aging 
of surfaces which is one of the most important rea- 
sons for pore formation. 

Depending on the electrochemical composition of Current 
Bursts the one or the other effect will dominate, and thus 
may lead to complete different surface morphologies when 
changing the electrochemical etching conditions; e.g. in 
a Si-HF-organic electrolyte system only by increasing the 
HF concentration (i.e. faster oxide dissolution and thus 
reduced influence of oxide) the electrochemical etching 
changes from electropolishing (strong oxide smoothing) 
over macropore formation (pores with diameters of sev- 
eral fim and smooth pore walls) to mesopore formation 



(strongly anisotropic, narrow pores with diameters of less 
than 400nm) [14]. 

Silicon has one of the most stable oxides of all semi- 
conductors. So for all other semiconductors the smoothing 
effect is reduced, leading generally to rougher surfaces and 
smaller pore diameters. In addition the surface aging of 
III-V compounds is more complicated, since there exist 
two different (111) surfaces; e.g. in GaAs only the {111} A 
planes (Ga-rich planes) appear as stopping planes. So in 
most III-V compounds not octahedra (eight (111) surfaces 
as stopping planes) as in silicon are etched but tetrahedra 
with only four (111) A surfaces as stopping planes, and 
the {111}B planes (As-rich) serve as preferential growing 
directions like the (100) directions in Si (see Fig. 7). 
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Figure 7: Tetrahedron-like pores oriented along < 111 > 
directions obtained in (lOO)-oriented n-GaAs (n = 10 17 
cm -3 ) at high current density (85 mA/cm 2 , galvanostatic) 
in HC1 (5%) electrolyte. 



7 Lateral Interaction of Pores 

In the case of stable pore growth active current bursts only 
exist at each pore tip, so there is no way that current burst 
can interact directly, e.g. by overlapping. Since several 
types of pores grow by forming chains of interconnected 
cavities, i.e. a self induced 

diameter modulation is an intrinsic feature of such 
pores, they can interact indirectly via the space charge 
region between pores. The serial resistance R of each pore 
can be described by 



R(t) 



A{t) 



(2) 



where p(t) is the effective specific resistance at a pore tip, 
l(t) the length of the pore and A(t) the chemically active 
area at the pore tip. Due to aging p(t) will periodically 
increase and decrease while the cavities are formed. The 
electric circuit consisting of a parallel connection of large 
numbers of individually oscillating pores is illustarted in 
Fig. 8. The consequences for the pore morphologies de- 
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Figure 8: A schematic representation of pores as oscillat- 
ing resistors in an equivalent cirquit diagram. 

pend strongly on the etching conditions as schematically 
illustrated in Fig. 9. At constant voltage condition (po- 
tentiostatic control) the current through each pore will be 
modulated; since the pores grow independently of each 
other there exists no phase coupling and the overall cur- 
rent is constant. 

For galvanostatic control we have to distinguish be- 
tween two cases: 

• At low pore densities (low external current density) 
there is no restiction for the diameter of the pores; 
i.e. the area A(t) may increase freely. 

• At high pore densities (high external current den- 
sity) the increase of the diameter is resticted by the 
neighboring pores; thus the area Ait) is limited. 

In the first case we will find diameter modulation (cf. Fig. 
9b)), but still no phase coupling between pores. In the 
second case the area of the pores can cot increase freely 
to compensate for the reduced effective specific resistance 
p{t). Thus the external voltage must increase to guarantee 
a constant external current. The voltage increase should 
synchronize the phases of growth for all pores. This is ex- 
actly what is found when, e.g., etching InP at very high 
current densities. In Fig. 10 a) self induced voltage os- 
cillations are shown which coincide with a simultaneous 



diameter increase of all pores (cf. the cross section of the 
pores in Fig. 10 b)). Obviously not all cycles of the pore 
growth have been synchronized globally allthough nearly a 
closed packed ordering of the pores was reached (cf. plane 
view of pore array in Fig. 10 c)). 
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Figure 9: A schematic representation of pore interaction 
explaining the morphologies observed in Si and III-V com- 
pounds; a) p oscillates and A is constant; b) p is constant 
and A oscillates; c) p and A oscillate simultaneously. Note 
that for galvanostatic conditions a constant total current 
is only possible for the case of random oscillator phases. 



8 Open-Loop-Control of Dynami- 
cal Systems 
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Figure 10: Data taken from an (100) oriented n-InP sam- 
ple with ni = 1.5 x 10 16 cm -3 anodized at a constant 
current density j = 100 mA/cm 2 . The electrolyte is again 
HC1 (5%). a) Voltage oscillations; b) cross-sectional SEM 
of the sample. The inset is the magnification of the nodes 
at the bottom of the porous layer; c) Induced by the next 
neighbor interaction, pores arrange into a close packed 
hexagonal 2D array. Correlation length is about 7 lattice 
constants. 



Contrary to the methods of Ott, Grebogi and Yorke [15] 
and of Pyragas [16], which allow for the stabilization of 
(if uncontrolled) unstable periodic orbits by parametric 
feedback, it may be possible to achieve stabilization even 
without feedback, i. e. by open-loop control. Applica- 
tion of non-feedback control to dynamical systems can be 
roughly classified into at least three classes: Addition of 
noise [17, 18] or of chaotic signals [19] may lead to suppres- 
sion of chaos. Second, a constant (but usually large) shift 
in a parameter may be utilized to shift a fixed point into 
a parameter regime where it is stable [20]. Third, peri- 
odic parametric perturbations can be used. Counterintu- 
itively, this is possible even by nonresonant perturbations 
[21]. The more straightforward approach is by resonant 
parametric perturbations [22, 23, 24, 25, 26, 27, 28] which 
in the simplest case is a weak sinusoidal signal, but can 
be any periodic signal. Which orbits can be stabilized 
by which signals in general can be investigated only if an 
explicit model of the system is given. 

9 Open-Loop-Control of Pore For- 
mation 

For technical applications, it would be desirable to gener- 
ate pores with a diameter of less than 500nm. However, 
at all etching conditions investigated up to now, in this 
diameter range always side-branching pores are obtained 
(cf. first row in Fig. 11). 

Fortunately, the experimental setup allows for an open- 
loop stabilization e.g. by modulation of the intensity of 
illumination or by modulation of the etching current. If 
the modulation frequency meets the intrinsic time scale 
of the pore formation (0.5 min modulation in row 1), the 
side pore formation is hindered by the reduction of the 
etching current at the right phase. Which is the relevant 
time scale (oxide dissolution time, passivation time, ...) 
depends on the etching conditions and the type of pores. 

Corresponding to the modulated current, the pores are 
regular, but with a modulated diameter (Fig. 11). The 
diameter modulation profile can be influenced to some ex- 
tent by the waveform, nota bene these waveforms are not 
identical. 

It should be noted, although the controlled state pos- 
sibly could be considered as an instable periodic orbit of 
the whole etching setup, the described control method yet 
does not use a feedback control scheme as the OGY scheme 
[15]. This would be a promising extension, however, it 
is expected that for the case of dendritic-like sidebranch- 
ing, the dynamics possibly cannot be reduced to a low- 
dimensional attractor. 




Figure 11: Branching of pores can be suppressed by an open- loop control method triggering to the system- inherent 
time scales (here: (lOO)-oriented n-Si, l-10^cm). The upper row (DMF, 4wt-% HF) corresponds to the case with less 
oxide generation. From middle row (H2O, 7wt-% HF) to lower row (H2O, 4wt-% HF), there is less oxide dissolution. In 
the first column the backside illumination current is constant; the triggering time decreases from 5 min in the second 
colume to 0.5 min in the last column. Here potentiostatic conditions (4V) are used; a 20% sinusoidal modulation 
(with the respective period) of the current density around its average value of 4 mA/cm 2 is realized by control of the 
amplitude of backside illumination. 



10 Conclusions 

The Current Burst model, a local stochastic and highly 
nonlinear model, and the Aging concept describing the 
time- and orientation-dependence of the passivation be- 
havior, provide on an intermediate level of abstraction 
a general approach to explain pore geometries, oscilla- 
tions and synchronization with a minimum of material- 
dependent, but experimentally accessible, parameters. The 
lateral interaction of Current Bursts gives rise to synchro- 
nization phenomena, and a percolation transition to global 
ordering. In a technologically relevant diameter range, 
where one otherwise experiences excessive sidebranching 
of pores giving dendritic-like structures, an open-loop con- 
trol can be successfully applied to suppress sidebranching 
of pores. To clarify quantitatively the dynamical processes 
in this system, and to gain detailed insight about the gen- 
eral controllability conditions for mesopores and perspec- 
tives for more sophisticated control strategies further ex- 
perimental and theoretical efforts have to be made. 
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